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ABSTRACT 

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological 
simulations of the ACDM cosmology and compare it to the shape of dark matter (DM) distribution and the 
shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation 
and stellar feedback (CSF) intracluster gas outside the cluster core (r > O.lrsoo) is more spherical compared to 
non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly 
oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses 
thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas 
cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected 
in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the 
CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in 
non-radiative simulations the difference between gas and potential shape at small radii is due to random gas 
motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock 
Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas 
are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters 
against observations can therefore constrain the amount of cooling in the intracluster medium and the presence 
of random gas motions in cluster cores. 



Subject headings: cosmology: theory, 
numerical 



1. INTRODUCTION 



galaxies: clusters: general - X-rays: galaxies: clusters - methods: 



In the prevailing, hierarchical cold dark matter (CDM) 
paradigm of cosmological structure formation, galaxy- and 
cluster-sized CDM halos are formed via accretion and merg- 
ing with smaller halos. The CDM paradigm predicts that DM 
halos are generally triaxial and are elongated along the direc- 
tion of their most recent major mergers. The triaxiality of DM 
halos has been demonstra ted in a nu mber of studies using nu- 
merical simulations (Frenk et al.ll988; Dubinski & Carlberg 
T9911 [Warren et al.||1 992; Thomas et al. 1998, Jing & Suto 
M5?'Suwa et al.'2003 ; Hopkins et al. 2005 ; Kasun & Evrar3l 
5D05 ; Allgood et al. 2006 ; Bett et al. 2007 ; Gottlober & Yepes 
^07 [ |Paz et al.|2008 ) and arises due to anisotropic accretion 
and mierging along filamentary structures. The degree of tri- 
axiality strongly correlates with the halo formation time (e.g., 
lAUgood et al.j2006l|Ho et al |2006[ |Wray et al.|2006, ), which 
impUes that at a given epoch more massive halos are more 
triaxial. For the same reason, triaxiality is sensitive to the lin- 
ear structure growth function and is higher in cosmological 



models in which halos form more recently (e.g., |Macci6 et al. 
[2008) . 

Although shapes of DM halos have been studied exten- 
sively in dissipationless A^-body cosmological simulations, 
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DM shape is difficult to probe observationally, though some 
handle on shape is provided by lensing studies ( e.g., Hoekstra 
let al.'2004','Parker et al.'2007:' Rozo etal.|2007[ [Evans & BrT 
die 2009 i Hawken & Bridle 2009]). Moreover, it is well known 
that including baryons in simulations modifies the shapes of 
DM halos, especially in th e case of significant gas dissi pation 
during galaxy formation ([Katz & G unn 1991 ; Evrard et al. 



1994 
MR 



Dubinski' 1 994', 'Kazantzidis et al. 2004; Springel et al. 
Ilayashi et al. 2007; Tissera et al. 2010, see Debattista 



et al. 2D08| and |Valluri et al. l| 2010| for discussion of the phys- 
ical nature of this effect). It is therefore of paramount im- 
portance to examine predictions for halo shapes using cosmo- 
logical simulations that include gas dynamics and dissipative 
processes accompanying galaxy formation. Further, the shape 
of the gas itself can be examined in such simulations and 
compared to the shape of the underlying potential, sourced 
predominantly by DM. Gas is expected to follow isopotential 
surfaces in hydrostatic equilibrium, so simulations may test 
whether cluster gas is in equilibrium on average and whether 
it can be used as a reliable tracer of the shape of the under- 
lying potential (B uote & Tsai 1995 ; Lee & Su to 2003 ; Flores] 
[et ar|2007[[Kawahara|2010p . 

Probing the shape of the gravitational pote ntial via gas, as 
was first suggested by Binney & Strimpel (1978| and ob- 
servatio nally tested by [Fabricant et al.,(1984j and [Buote &| 
[Caniza res ( 1996), can open mteresting avenues for using the 
shapes of DM halos around observed galaxy clusters to both 
test the CDM paradigm and constrain the amount of halo gas 
that dissipated and was converted into stars during halo for- 
mation. This is particularly relevant for galaxy clusters, where 
high-quality X-ray imaging data now exists for large samples 



of clusters. |Kawahara[ ( 20 1 ) recently analyzed axis ratios of 
X-ray clusters from the XMM-Newton catalogue and found 
relatively good agreement with the CDM predictions of [Jing] 
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& Suto ( 2002 1 based on dissipationless simulations of a large 
clus ter sample, confir ming findings of Buote & Tsai ( 1995]! ^ 
and'Flores et al. based on a single simulated systemFJ 

[Buote & Tsai ( 1995 ) were also the first to show that for the 
simulated cluster they studied the shape of the X-ray isophotes 
reflected the shape of the underlying three-dimensional gas 
distribution and potential. 



Table 1 

Properties of the Simulated Clusters at z = 



Fang et al. ( 2009) compared the ellipticities of X-ray sur- 
face brightness isophotes of clusters simulated with and with- 
out radiative cooling and star formation. Focusing on a single 
cluster from the sample we analyze in this paper, they showed 
that the shape of gas can be quite flattened when gas cools 
significantly and settles into rotating thick disk. They also 
showed tha t this flattening is d etectable in the shape of X-ray 
isophotes. Fang et al. (20091 also argued that the flattened 
shape of the gas distribution m the simulated cluster implies 
that gas does not trace potential in the inner regions due to 
rotational support. Their results therefore demonstrate that 
shapes of X-ray isophotes in cluster cores are a useful diag- 
nostic of amount o f cooling and gas motions in cluster cores. 
Fang et al. ( 2009[ ) have also compared isophote shapes for 
synthetic Chandra observations of a sample of clusters simu- 
lated with cooling to observations and concluded that elliptic- 
ity profiles in runs with radiative cooling do not match obser- 
vations. They attributed the discrepancy to ongoing, signifi- 
cant cooling in the cores of simul ated clusters that is absen t 
from the cores of real clusters (e.g., [Peterson & Fabian |2006l l. 

In this paper we present analysis of the three-dimensional 
shapes of intracluster gas, DM, and underlying gravitational 
potential using high-resolution cosmological simulations of 
galaxy clusters formed in the ACDM cosmology. Our work 



extends the work of [Fang et aL ( 2009 1 by presenting more 
detailed analysis of three-dimensional shape profiles of gas, 
DM, and gravitational potential for the full sample of clus- 
ters. In addition, we focus on the effects of cumulative cool- 
ing and dissipation during the entire cluster evolution on the 
shape of potential and gas distribution at intermediate radii 
(r > 0.2r5 oo), where dissipation ma kes potential more spher- 
ical (e.g., Kazantzidis et al."2004'). This effect was not in- 
vestigated in Fang et al. (2009). We show explicitly that 
dissipation leads to more spherical shapes for ICM gas out- 
side cluster cores (0.1 < r/rsoo ^ reflecting the corre- 
sponding effect on the DM distribution. We also find that the 
shape of gas matches the shape of the gravitational potential at 
these radii in general, but deviates from it at smaller radii and 
at r > rsoo where assumption of the hydrostatic equilibrium 
breaks down. At smaller radii (r < O.lrsoo) gas distributions 
in simulations with cooling become oblate, reflecting the par- 
tially rotation-supported thick disks into which gas settles as 
it loses its thermal energ y by cooUng, a re sult that is qual- 
itatively consistent with Fang et al. (2009 1. We predict gas 
shapes as may be determined observationally by estimating 
ellipticities from mock X-ray maps of the same clusters and 
show that one can constrain cluster gas physics by comparing 
ellipticity profiles of simulations with and without dissipation 
to those of observations. 

This paper is organized as follows. In Section [2| we de- 
scribe our cluster simulations. In Section [3] we describe the 



^ A similar test of the shape determined using t he observed project ed 
galaxy distribution of clusters have been presented by |Plionis et al.|pQ06| . 

* Here and throughout this paper rsoo denotes the cluster-centric radius 
enclosing a mean overdensity of 500pc(z), where pdz) is the critical density 
of the universe at the redshift of analysis. 
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method of estimating axis ratios and the results for the three 
dimensional shapes of clusters. We give our estimates for el- 
lipticities of mock X-ray maps as an observable prediction in 
Section [4[ We provide a summary and discussion of our re- 
sults in Section [5] 

2. THE SIMULATIONS 

We analyze high-resolution cosmological simulations of 16 
cluster-sized systems in a flat ACDM model: ilm = 1 - iljy = 
0.3, fib = 0.04286, h = 0.7 and erg = 0.9, where the Hubble 
constant is defined as 100/; km s~' Mpc~\ and is the mass 
variance within spheres of radius 8/i~' Mpc and serves to nor- 
malize the power spectrum. The simulations were performed 
using the Adaptive Refinement Tree A^-b ody+gasdynamics 
code ( |Kravtsov|1999t Kravtsov et al.|2002] l, an Eulerian code 
that uses adaptive refinement in space and time, and (non- 
adaptive) refinement in mass (Klypin et al. 2001 ) to reach the 
high dynamic ranges required to resolve the cores of halos 
formed in self-consistent cosmological simulatio ns. The sim- 
ulations pres ented here are disc ussed in detail in [Nagai et al. 
( 2007a[ l and Nagai et al. (2007b| and we refer the reader to 
these papers for more details. Here we summarize the rele- 
vant parameters of the simulations. 



In order to assess the effects of gas cooling and star for- 
mation on the cluster shapes, we conducted each cluster sim- 
ulation with two different prescriptions for gasdynamics. In 
one set of runs we treated only the standard gasdynamics for 
the baryonic component without either radiative cooling or 
star formation. We refer to these as non-radiative (NR) runs. 
In the second set of runs, we included gas cooling and star 
formation (CSF). In the CSF runs, several physical processes 
critical to various aspects of galaxy formation are included: 
star formation, metal enrichment and thermal feedback due to 
supernovae Type II and Type la, self-consistent advection of 
metals, metallicity-dependent radiative cooling and UV heat- 
ing due to a cosmological ionizing background (see Nagai] 
|et al.|[2007b[ for details of the metallicity-dependent radia- 
tive cooling and star formation). These simulations there- 
fore follow the formation of galaxy clusters starting from 
well-defined cosmological initial conditions and capture the 
dynamics and properties of the intracluster medium (ICM) 
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in a realistic cosmological context. However, some poten- 
tially relevant physical processes, such as active galactic nu- 
clei (AGNs) bubbles, magnetic fields, and cosmic rays, are 
not included. Consequently, the simulated cluster properties 
have limited application to real systems, most notably in the 
innermost cluster regions, where these processes are likely to 
be more important. 

We ran our high-resolution simulations using a uniform 
128^ grid and eight levels of mesh refinement in computa- 
tional boxes of 120/1"' Mpc comoving on a side for CLIOI- 
107 and 80 /i"' Mpc for CL3-24. These simulations achieve 
a dynamic range of 32768 and a formal peak resolution of « 
3.66/1"' kpc and 2.44 /1"' kpc, corresponding to an actual reso- 
lution of w 7 /1"' kpc and 5 /1"' kpc, for the 120 and 8O/1"' Mpc 
boxes, respectively. Only regions of ^ 3 - 10 /1"' Mpc sur- 
rounding each cluster were adaptively refined. The remain- 
ing volume was followed on the uniform 128^ grid. The 



particle mass. 



corresponds to an effective 512 parti- 



cles in the entire box, or a Nyquist wavelength of Any = 
0.469/!"' Mpc and 0.312/i"' Mpc comoving for CLlOl-107 
and CL3-24, respectively. These correspond to 0.018 /?"' Mpc 
and 0.006/?"' Mpc in physical units at the initial redshifts of 
the simulations. The DM particle mass in the region around 
each cluster was mp ~ 9.1 x lO'*/i"'M0 for CL 101-1 07 and 

nip ~ 2.7 X 10^/i"'Mo for CL3-24, while other regions were 
simulated with lower mass resolution. For this paper, we only 
report values at cluster-centric distances larger than 0.03r50() 
where all clusters are well resolved. 

In order to test our simulation results against observations, 
we created mock Chandra X-ray images along three orthog- 
onal projections for each simulated cluster To minimize sta- 
tistical fluctuation due to the Poisson noise, each image has 
an exposure time of 100 ks, corresponding to deep X-ray ob- 
servations. Instrumental responses of Chandra were included 
in the mock image data. An overview of the methods used to 
generate the mock images is given in Section [4] Detailed de- 
scriptions of these mock images can be foundin Section 3.1 
ofNagaietal. (2007b). 

To investigate the dependence of gas and halo shapes on 
the dynamical states of the clusters in our simulation set, we 
divided our sample into relaxed and unrelaxed clusters based 
on visual examinations of their mock X-ray images. A typical 
relaxed cluster is an object in which all three of its orthogonal 
images have a regular morphology and prominent substruc- 
tures are absent. Details of the classification can be found in 



Nagai et al |(l2007a[ ). 

In Table llj we report M500 (the mass within rsoo), r^oo, and 
our classification of each of the X-ray images along the three 
orthogonal projections as relaxed or unrelaxed for our sample 
of 16 z = clusters. 

3. THREE DIMENSIONAL SHAPES 
3.1. Methods 



We esti mate axis ratios fo l lowin g Dubinski & Carlberg 
(|I99I[) and [Kazantzidis et al.|(2004[). For DM particles, we 



compute the mertia tensor 



(1) 



where Xj is the / coordinate of particle a, Wq is the particle 
mass, and the sum is over all particles in a shell of width Ar. 
The principal axis lengths are obtained by diagonalizing /,j 
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Figure 1. Average ellipsoidal axis ratio profiles for the relaxed z = clusters 
from the CSF run (left panels) and the NR run (right panels). The upper 
panels show the profiles for the short-to-long axis ratio c/a, and the bottom 
panels show the profiles for the intermediate axis ratio b/a. In all panels, the 
solid line corresponds to dark matter (DM), the dashed line corresponds to 
gas, and the dot-dashed line coiTesponds to gravitational potential. The error 
bars show Icr error on the mean axis ratio for gas. The magnitude of the 
errors on the mean axis ratios for gas is similar to those of DM and potential. 



using an iterative scheme. We begin by computing the in- 
ertia tensor in spherical shells and computing the axis ratios 
q = b/a and s = c/a <q using the eigenvalues of the inertia 
tensor (we adopt a > /? > c by convention). In the subse- 
quent iterations, we compute principal axes at a given radius 
r by summing over particles within ellipsoidal shells of width 
Ar, defined using orientations and values of the principal axes 
from the previous iteration. The generaUzed elUpsoidal dis- 
tance of particle a is 




(2) 



where the primed coordinates are particle positions rotated 
into the frame of the principal axes of the inertia tensor. The 
centers of all ellipsoids are fixed to be the center of the cluster, 
defined as the location of the most bound particle. This pro- 
cess is continued until convergence of q and s to better than 
1%. 

Note that Dubinski & Carlberg (1991) weighted each term 
in their calculation of the inertia tensor by r^/^, to mitigate 
the influence of substructures, prevalent at the outskirts of ha- 
los, on axis ratios comp uted within a certain radius. As in 
Kazantzidis et al. ( 2004| ), we find that this weighting makes 
only a small difference for axis ratios computed within nar- 
row radial bins so we do not employ this weighting in our 
analysis. In this paper we use the mean r^ (Equation (|2]))of 
particles within each ellipsoidal shell (equivalent to the length 
of the major axis of the shell) as our measure of cluster-centric 
distance, unless stated otherwise. 

We estimate the axis ratios for gas in a similar way with 
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Figure 2. Difference in tlie axis ratio profiles of dark matter (left panel), gas (middle panel), and gravitational potential (right panel) between CSF and NR runs 
averaged over relaxed clusters at = 0. The difference is defined as the CSF axis ratios minus the NR axis ratios. The black solid line is the difference for the 
short-to-long axis ratio, (A(c/a)). The red dashed line is for the middle-to-long axis ratio, (A(i/o)). 



weights Wa = Pgas,aKi> where Pgas.a and are the gas density 
and volume of the a-th grid cell. We estimate the axis ratios of 
the surfaces of constant gravitational potential by computing 
the inertia tensor of all cells with potential within some range 
[<i>, $ + A(l)], taking = 1 for all such cells. 

Large subhalos can bias the axis ratios in a given radial 
bin to lower values, generating a local fluctuation in the axis 
ratio profile. To minimize fluctuations due to substructures 
we remove particles bound to subhalos of masses Msub > 
1O'^/i~'M0. The identification of subhalos and bound parti- 
cles follows the procedure described in Kravtsov et al., ( ,2004^ . 



3.2. Results 

Figure [T| shows the axis ratio profiles for dark matter, gas, 
and gravitational potential averaged over the relaxed clusters 
at z = for the CSF and NR runs respectively. Results are 
similar for unrelaxed clusters but with considerable scatter. In 
both CSF and NR runs, gravitational potential is much more 
spherical than DM. The potential of a given thin triaxial mass 
shell with axis ratios q and s (the homeoid) is constant within 
the shell and has a triaxial shape outside of isopotential sur - 
faces defined by the ellipsoid (e.g., Binney & Tremaine 2008 1: 



r 



T+l T + q^ 



T + S^ 



where r is the label of the surface. The shape of the isopoten- 
tial surfaces corresponding to a given ellipsoidal shell there- 
fore becomes more spherical as the distance from the shell 
increases. Total potential at a given distance from the cluster 
is the superposition of the potentials generated by ellipsoidal 
shells within this distance and a constant potential generated 
by ellipsoidal shells outside. The shape of the isopotential sur- 
faces therefore will always be more spherical than the shape 
of the underlying mass distribution that gives rise to the po- 
tential. 

Figure [1] shows that at r > O.lrsoo the shapes of DM, grav- 
itational potential, and gas are more spherical in the relaxed 
CSF clusters than in the relaxed NR clusters, an effect iden- 
tified previously in several studies (see Section [TJ. This ef- 
fect appears to be due to the fact that potential becomes more 
spherical at large radii for a more concentrated mass distribu- 
tion resulting from baryon dissipation and adiabatic response 



of the parti cle orbits to such c hange of potential ( [Debattista] 
et al. 2008 ; Valluri et al.|2010| l. 

The effect is more apparent in Figure l2] which shows the 
differences in axis ratios between the relaxed CSF and NR 
clusters defined as 



A(c/fl)= (c/a)p5p-(c/ 



(3) 



and similarly for b/a. It is evident that baryonic dissipation 
causes relaxed DM halos to become significantly more spher- 
ical in their inner regions, an effect that remains significant 
out to r5oo and beyond. The average axis ratio shifts drop 
from (A(c/fl)) ~ {A(b/a)) - 0.3 at O.lrjoo to {A(c/a)) - 
(A(c/a)) ^0.1 at r^oQ. Changes in c/a and b/a are very 
nearly the same at all radii. 

For gas, the average change in both axis ratios is 
(A(c/fl)) w {A(b/a)) ss 0.1 at 0.2r5()o, and is decreasing 
slowly to zero at ^500. However, at smaller radii, r < O.Irsoo, 
there is a positive change in the intermediate axis ratio b/a 
and a negative change in the short-to-long axis ratio c/a. At 
r w O.OSrsoo, {A{c/a)) w -0.3 and continues to decrease with 
radius, while {A{b/a)) increases to w 0.2 at 0.06r5oo and de- 
creases inward. The significant decrease in c/a indicates that 
the gas assumes an oblate shape in the inner regions of the re- 
laxed CSF clusters compared to the prolate gas shapes in their 
NR counterparts. This is consistent with the expectation that 
ongoing cooling in the inner cluster region leads to formation 
of a thick oblate disk (Fang et al. 2009 ). 

The difference between the shapes of the gravitational po- 
tential of relaxed CSF and NR clusters is qualitatively sim- 
ilar to that of the DM. The average change in axis ratios is 
(A(c/fl)) w {A{b/a)) w 0.2 at r O.lrsoo, and decreases to 
nearly zero at rsoo. There is little difference between {A(c/a)) 
and {Aib/a)). 

Finally, Figure [T] shows that gas traces the shape of gravi- 
tational potential ror relaxed NR clusters in the radial range 
0.06 < r/rsoo < 1, indicating that gas is in approximate hy- 
drostatic equilibrium within the potential. However, in the 
relaxed CSF clusters the gas and potential shapes only match 
at over a relatively narrow range of radii, 0.2 < r/rsoo < 0.4, 
signaling departures from hyd rostatic equilibrium. We dis- 
cuss this further in Section [3T3] 

Figure [3] shows c/a versus b/a at r/rsm = (0.1, 1.0,2.0) for 
our relaxed and unrelaxed cluster samples. The average c/a 
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Figure 3. Minor-to-major axis ratio (cjd) vs. intermediate-to-major axis ratio (bja) at r = O.lrjoo (l^ft panels), r = rsoo (center panels) and r — 2r5oo (right 
panels) for DM (top panels), gas (middle panels), and gravitational potential (bottom panels) at z = 0. Solid points are relaxed clusters and open points are 
unrelaxed clusters. CSF clusters are represented by circles, NR clusters are represented by squares. The dashed line is c = i. 



values are summarized in Table |2] As could be expected, the 
DM distribution in unrelaxed clusters is more triaxial on av- 
erage compared to relaxed clusters. The DM halos in the NR 
clusters become more triaxial at smaller cluster-centric radii, 
consistent with previous findings in dissipationless simula- 
tions (e.g., Allgood et al.|2006 ). Conversely, the CSF clusters 
are rounder at small radii due to the effects of baryonic dis- 
sipation on halo shapes which are m ost prominent near halo 
centers (e.g. Kazantzidis et al.|2004| l. 

Intracluster gas is very spherical at large radii in almost all 
relaxed NR and CSF clusters. Unrelaxed clusters tend to have 
more triaxial gas distribution compared to relaxed clusters. A 
similar trend can be seen for potential. 

Figure [4] shows cja as a function of cluster mass Msoo- 
There is at most only a weak trend of decreasing cja with 
cluster mass for both the CSF and NR relaxed clusters seen in 



DM and potential. This is consistent with the results of stud 
ies based on larg e statistical samples of halos in dissipation 
less simulations (Kasun & Evrardl2005' Allgood et al."2006 
Gottlober & Yepes 2007; Ragone-Figueroa & Plionis 2007 
Maccio et al. 2008 ) which find cja cx M-'""-'-" "^J. Our sam- 
ple of clusters is too small to detect such a trend. 

3.3. Comparing the Shape of Intracluster Gas to the Shape of 
Gravitational Potential 

In hydrostatic equilibrium, the isodensity surfaces of gas 
should trace the isopotential surfaces. Consequently, any 
differences between the measured shapes of gas and poten- 
tial indicate deviations from hydrostatic equilibrium. Fig- 
ure [T] shows that such differences are present at small radii 
r < 0.06r5oo in relaxed NR and r < 0.2r5()o in relaxed CSF 
clusters, and at larger radii r > r^oo and > 0.6r5oo in the NR 
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Figure 4. Minor-to-major axis ratio c/ci evaluated at r = O.lrsoo (left panels), rjoo (center panels), and 2rjoo (right panels) as a function of Mjoq. DM shape 
profiles are shown in the top row of panels, gas in the middle row, and gravitational potential in the bottom row. All cluster profiles are computed at z = 0. Solid 
points are relaxed clusters, and open points are unrelaxed clusters. CSF clusters are represented by circles, NR clusters are represented by squares. 
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Figure 5. Relative contributions of random and rotational motions support- 
ing gas against gravity in the relaxed CSF and NR clusters. The left panel 
shows the ratio of the gas velocity dispersion to the circular velocity as a 
function of cluster-centric position. The right panel shows the ratio of the 
rotational velocity to the circular velocity. 



and CSF runs, respectively. All of these differences are due 
to the presence of gas bulk motions, although the nature and 
origin of these motions is differ ent at small and large radii and 
in NR and CSF runs (see |Laue t al. 2009 ). 

Figure [5] shows the ratio of the isotropic gas velocity dis- 
persion CTgas and gas rotational velocity (calculated within 
each ellipsoidal gas shell) to the circular velocity (defined as 
Vciic = \/ GM(< r)/r, where r is the mean ellipsoidal radius 
of the shell (Equation Q). This ratio indicates the relative 
contribution of random gas motions in support against grav- 
ity. For gas in perfect hydrostatic equilibrium CTgas should be 
zero. However, we see that in both the NR and CSF clus- 
ters CTgas ~ 0.2-0.4vcirc for 0.1 < r/rsQo ^ 1, corresponding 
to a fraction of pressure support due to ranclom gas motion of 
about 16%. These motions are dynamical in origin (due to 
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Table 2 

axis ratio c/a of the 16 simulated clusters 







z = 











0.6 






0.1 


0.5 


1.0 


2.0 


0.1 


0.5 


1.0 


2.0 


All A Jf 

All DM 


CSF 


0.66 it 0.05 


0.47 ±0.05 


0.54 ± 0.04 


0.39 ±0.06 


0.59 ±0.06 


0.32 ±0.06 


0.33 ±0.07 


0.36 ± 0.06 




NR 


0.37 it 0.03 


0.38 ±0.04 


0.41 ± 0.05 


0.28 ±0.05 


0.40 ± 0.04 


0.38 ±0.04 


0.41 ± 0.04 


0. 14 ± 0.05 


potential 


CSF 


0.59 it 0.10 


0.58 ±0.08 


0.55 ± 0.06 


0.53 ±0.05 


0.59 ±0.10 


0.50 ±0.07 


0.50 ±0.07 


0.38 ±0.05 




NR 


0.49 it 0.08 


0.57 ±0.07 


0.56 ±0.05 


0.49 ±0.05 


0.33 ±0.08 


0.44 ±0.06 


0.54 ±0.06 


0.22 ±0.07 


gas 


CSF 


0.58 it 0.05 


0.67 ±0.67 


0.74 ±0.04 


0.70 ±0.06 


0.68 ± 0.04 


0.69 ±0.05 


0.69 ±0.06 


0.59 ±0.05 




NR 


0.58 it 0.05 


0.65 ±0.05 


0.67 ±0.06 


0.63 ± 0.06 


0.39 ±0.07 


0.55 ±0.06 


0.66 ±0.04 


0.26 ±0.08 


Relaxed DM 


CSF 


0.76 it 0.02 


0.63 ±0.05 


0.67 ±0.04 


0.46±0.10 


0.51 ±0.12 


0.39±0.10 


0.40 ±0.09 


0.39 ±0.1 3 




NR 


0.41 ±0.02 


0.46 ±0.06 


0.48 ±0.07 


0.30±0.10 


0.32 ±0.04 


0.42 ±0.03 


0.46 ±0.01 


0.17±0.11 


potential 


CSF 


0.91 ±0.01 


0.82 ±0.02 


0.71 ±0.02 


0.65 ±0.03 


0.82 ±0.02 


0.70 ±0.02 


0.69 ±0.03 


0.49 ±0.07 




NR 


0.65 ±0.03 


0.72 ±0.06 


0.67 ±0.05 


0.58 ±0.02 


0.36±0.12 


0.64 ±0.02 


0.69 ±0.04 


0.22 ±0.1 3 


gas 


CSF 


0.61 ±0.07 


0.79 ±0.05 


0.80 ±0.03 


0.82 ±0.02 


0.68 ±0.06 


0.78 ±0.02 


0.80 ±0.04 


0.75 ±0.05 




NR 


0.61 ±0.07 


0.74 ±0.03 


0.76 ±0.02 


0.72 ±0.05 


0.43 ±0.13 


0.67 ±0.04 


0.72 ±0.05 


0.28 ±0.1 8 



accretion of gas and motions of cluster galaxies) and are not 
affected by cooling. 

The gas motions do not affect the shape of the gas signifi- 
cantly at radii where (Xgas/vcii-c ^ 0.4, but have significant ef- 
fect for larger values of cTgas/vdrc, as is apparent from the dif- 
ferences between the shapes of gas and potential at small and 
large radii in the NR clusters. The net effect of these motions 
is relatively rounder gas distributions compared to the shapes 
of isopotential surfaces. 

In addition to random motions, the gas in the CSF runs 
exhibits significant ordered rotational motions at r < O.Srsoo- 
These motions arise due to the angular momentum of the ICM 
gas, which leads to rotation as gas cools and contracts. These 
motions are responsible for the deviations between the shapes 
of gas and potential at these radii. As we can see in Figure [T| 
the effect of the ordered motions is opposite to that of random 
motions, namely, rotational motions lead to gas distributions 
that are more flattened compared to isopotential surfaces. 

The effect of cooling, and the ordered rotational motions 
that result from cooling, manifests as a rapid decrease of c/a 
at small radii, even as b/a remains approximately constant. 
The effect of random gas motions results in a rapid increase 
of both the c/a and b/a ratios. Measurements of the ellip- 
ticity profiles in clusters can therefore constrain the cooling 
of gas and magnitude of residual gas motions in their cores. 
Although X-ray spectroscopy is a much more sensitive tool 
to constrain the contemporary cooling, rotational motions and 
their effect on gas ellipticity also constrain the net cooling 
that has been occurred in the past. They can therefore pro- 
vide potentially useful and complementary constraints on the 
thermal history of gas in cluster cores over the past several 
billion years. However, in order to make sure that such con- 
straints are feasible we must check that the trends observed in 
the three-dimensional gas distribution are evident in the X-ray 
photon maps of these clusters. We present such an analysis in 
the next section. 

4. ELLIPTICITY OF X-RAY IMAGES OF SIMULATED CLUSTERS 
4.1. Methods 

Although one can expect that at radii where gas traces grav- 
itational potential the X-ray isophotes can be used to quantify 
its shape, it is not immediately clear whether this can be done 
with the accuracy sufficient to detect the effects of cooling 
discussed in the previous section. In order to compare our re- 
sults on intracluster gas shapes in Section |3] to observations, 



we estimate the ellipticities of mock X-ray photon maps gen- 
erated from the same z = simulated clusters. Below we give 
a brief overview of the mock X-ray maps. We refer the reader 
to Section 3.1 of |Nagai et al.| ( |20Q7b| for a detailed descrip- 
tion. 

First we create X-ray flux maps of the simulated clusters 
viewed along three orthogonal projections. The flux map is 
computed by projecting the X-ray emission of hydrodynamic 
cells enclosed within 3rvii of a cluster along the line of sight. 
The X-ray emissivity in each computational grid cell is com- 
puted as a function of proton and electron densities, gas tem- 
perature and metal abundance. Emission from gas with tem- 
perature less than 10^ K is excluded as it is below the Chandra 
bandpass. We then convolve the emission spectrum with the 
response of the Chandra front-illuminated CCDs and draw 
a number of photons at each position and spectral channel 
from the corresponding Poisson distribution. Each map has 
an exposure time of 100 ks (typical for deep observations) 
and includes a background with the intensity corresponding 
to the quiescent backgro und level in the ACIS-I observations 
( [Markevitch et al.|2003| l. The resolution of all the maps is 6 
kpc pixer^ We use at least 25 pixels per bin for ellipticity 
measurements. 

From these data, we generate images in the 0.7-2 keV 
band and use them to identify and mask out all detectable 
small-scale clumps, as is routinely done in observational stud- 
ies. Our clump detection is fully automatic and based on the 
wavelet decomposition algorithm ( Vikhlinin et al. 1998). The 
holes left by masking out substructures in the photon map are 
filled in by the values from the decomposed map of the largest 
scale in wavelet analysis. We have tested that this method pre- 
serves the global shape of the photon distribution well. The 
background is removed when estimating ellipticities as it can 
bias them low at radii where background dominates the in- 
trinsic emission. Throughout this paper we assume the cluster 
redshift is Zobs = 0.06 for the z = sample. 

We define the ellipticity as 

b 

e=l--, (4) 
a 

where a and b are the semi-major and the semi-minor axes 
of the projected ellipse respectively. The ellipticities of the 
X-ray photon distributions are determined using the same al- 
gorithm, based on the inertia tensor, as was used for the three- 
dimensional shapes in Section [3] Instead of using particle 
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mass, we calculate the inertia tensor using weights given by 
the photon counts in each map pixel. We have estimated the 
ellipticity within both differential radial bins, and cumulative 
ellipticity within a given radius. In addition we have estimated 
ellipticity within radial shells defined by isophotes with differ- 
ent flux levels (isophotal bins). We find from visual compar- 
ison of the X-ray contours and the fitted ellipses that radial 
bins give reliable ellipsoidal fits to the X-ray isophotes (see 
Fig. [6]for an example). 

4.2. Test Case: CL7 

We test our method for estimating ellipticity using CL7, one 
of the most relaxed clusters in our sample. Figurel5|shows the 
mock X-ray maps for the y-projection of the NR and CSF runs 
(labeled as CL7:NR and CL7:CSF respectively). Also shown 
are the isophotal contour lines (red) and the fitted ellipses of 
the photon distribution (green) derived using the method de- 
scribed in Section l4n 

These images show that the ellipticities in the outer regions 
are similar for the NR and CSF runs, while they are very dif- 
ferent in the core (< O.lrgoo). The ellipticity increases to- 
wards small radii for the CSF run, while the opposite trend 
is seen for the NR run, where the isophotes are becoming 
slightly more spherical in the inner regions. These images 
also demonstrate that the best-fit ellipse (indicated with green 
contour) describes the actual photon distribution (indicated 
with red contours) well at all radii, including the flattened gas 
disk in the inner most regions of the CSF run. 

Figure |7] shows the ellipticity profiles of the CL7:NR and 
CL7:CSF runs viewed along three orthogonal projections. 
Note that we have scaled the radius to r = ^fah (where a and 
b are the semi-major and se mi-minor axes of the fitted el- 
lipse), to be consistent with Fang et aL] ( |2009 1. In the NR 
run, the ellipticity profiles for all three projections are quite 
similar over the range r > O.lrgoo, with nearly constant el- 
lipticities of e w 0.2 (though it decreases to e w 0.05 at rgoo 
in the jc-projection). For the CSF run, two projections have 
very similar ellipticity profiles (the y- and z-projections), with 
e w 0.2 at r w rgoo and increasing towards smaller radii, reach- 
ing e w 0.4-0.45 at r w O.lrgoo. The x-projection, on the 
other hand, exhibits low ellipticity comparable to the NR run, 
e « 0.1, and shows no strong trend with radius. These re- 
sults are consistent with the picture that the inner regions of 
the CSF run consist of a disk-like structure in the core, which 
appears highly elliptical when viewed edge-on (the y- and z- 
projections), while the flattened disk appears more spherical 
when viewed face-on (in the x-projection). 

Recently, the same set of simulated clusters are analyzed 



dFang et al.|2009| hereafter F09). In particular, the ellipticity 



prohles of CL7 was presented in detail (cf. Figure 7 in F09). 
Although our results are in qualitative agreement with those of 
F09 as to the formation of flattened gas structure in the cluster 
core, we find substantial quantitative differences regarding the 
impact of the baryonic dissipation. For example, we find that 
the ellipticity profiles for the CL7:CSF run are generally lower 
than the corresponding F09 profiles. In the y-projection at 
r = 0.3r5oo, we find e « 0.2 as opposed to the F09 value of 
e « 0.5. Our ellipticities are even significantly (by « 0.2- 
0.3) smaller in the cluster core. In addition. Figure [7] shows 
that significant flattening of the isophotes due to cooling is 
confined to r < 0.2r5oo, not out to 0.4r5oo as stated by F09. 

We have also checked ellipticities from cumulative bins 
rather than annular bins and found that the ellipticity increases 



only slightly to e = 0.3 at r = 0.3r5oo. By inspection, the 
F09 surface brightness map for the y-projection of CL7:CSF, 
shown in their Figure 1, appears to be inconsistent with e = 0.5 
at r = 0.3r5oo, but is consistent with the isophotes and our fits 
shown in Figures [6| and [7] To pin down the discrepancy, we 
focus on the y-projection of CL7:CSF and use our elliptic- 
ity code on both our FITS file and the FITS file given by the 
authors of F09. For consistency, we use cumulative bins in 
both cases. We find that the ellipticity profiles derived from 
the two different FITS files have very similar shape. The only 
difference is that the ellipticity profile from their FITS file ap- 
pears to be shifted to larger radii by a factor of 1.75. Upon 
inspection of the two FITS file, we find that the length scale 
of their FITS file is 1 .75 larger compared to our own FITS file. 
The origin of this 1.75 factor is unexplained in F09. We fur- 
ther find that the ellipticity profile estimated using our code 
on their FITS file is nearly identical to the profile for the y- 
projection shown in Figure 7 of F09. Thus we conclude that 
our discrepancy with F09 is most likely due to this 1.75 factor 
in their FITS files. 

Therefore, although we see effects qualitatively similar to 
those pointed out by F09, the actual magnitude of the effect 
of dissipation of ellipticity of the ICM gas is much smaller 
than measured by F09 and is confined to the inner r < 0.2r5oo 
of the clusters. 

4.3. Results 

Figure |8] shows the ellipticity profiles derived from mock 
X-ray maps averaged over the three orthogonal projections 
separately for all clusters as well as subsets of the images of 
relaxed clusters and the images of unrelaxed clusters (see Ta- 
ble [T]). The radial coordinate here is actually the semi-major 
axis a of the fitted ellipse in units of r^Qo. We note explicitly 
that the three-dimensional results in Section[3]were quoted as 
axis ratios, while the results of this section are given in ellip- 
ticities (e = 1 -b/a). 

The Figure shows that X-ray isophotes are more flattened 
in cluster cores in the CSF runs compared to the NR clusters. 
There is a clear rapid upturn in ecsF at r < O.lrgoo reflecting 
rotational motions of gas in these runs. Note that we do not 
confirm results of F09 who claimed significant flattening of 
isophotes due to rotation out to r « 0.4r5oo based on the anal- 
ysis of the same simulations we use in this study. 

There is a downturn in e for the NR clusters at similar 
radii (r < O.lrgoo) reflecting the effects of random gas mo- 
tions. These trends are consistent with the results for the 
three-dimensional distributions presented in Sectionl3]and are 
just as pronounced. The difference between CSF anoNR runs 
is particularly large for relaxed clusters, which have had more 
time since their most recent major mergers during which gas 
cooling could proceed. 

Outside the cluster cores, 0.1 < r/rgoo < 1, the ellipticities 
are approximately constant in both the CSF and NR runs. The 
average ellipticities in this range differ for the two sets of re- 
laxed clusters. Relaxed clusters in the CSF simulations yield 
an average ellipticity of (ecsf) ~ 0.1 -0.15 while the relaxed 
NR clusters have (eNk) ~ 0.15-0.25. The more spherical DM 
distribution and shape of the isopotential contours in the CSF 
clusters compared to the NR clusters are clearly discernible 
in X-ray maps of the relaxed systems. Note that these differ- 
ences are highly significant both in individual bins for relaxed 
clusters, and because they persist over a wide range of radii. 
The average ellipticities of unrelaxed clusters in the CSF and 
NR clusters are (ecsf) ~ 0.25 and (cnr) ~ 0.25-0.35 respec- 




Table 3 

Ellipticity of the 48 mock x-ray maps for the z = clusters 



r/rsoo = 


0.05 


0.3 


1.0 


2.0 


All 


CSF 


0.27 ±0.03 


0.20 ±0.02 


0.17±0.02 


0.40 ±0.04 




NR 


0.19 ±0.03 


0.30 ±0.03 


0.20 ±0.02 


0.38 ±0.03 


Relaxed 


CSF 


0.32 ±0.06 


0.12±0.01 


0.14±0.02 


0.32 ±0.06 




NR 


0.17 ±0.05 


0.21 ±0.02 


0.13±0.01 


0.31 ±0.04 


Unrelaxed 


CSF 


0.25 ±0.03 


0.27 ±0.04 


0.20 ±0.03 


0.43 ±0.05 




NR 


0.21 ±0.04 


0.37 ±0.06 


0.25 ±0.03 


0.45 ± 0.04 
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Figure 7. Ellipticity profiles of the CL7 cluster with NR (left panel, 
CL7:NR) and in the CSF run (right panel, CL7:CSF) viewed along three 
orthogonal projections (.v, y, z)- 



lively. The mean cluster ellipticities and their la errors on 
the mean are summarized in Table |3] Note that the scatter is 
substantial and constraining the net effect of cooling and the 
prevalence of random motions will require averaging over a 
sample of relaxed clusters. 

5. SUMMARY AND DISCUSSION 



We have investigated the axis ratios of DM mass, hot gas, 
and gravitational potential using 16 clusters simulated from 
the same initial conditions, but with different baryonic physics 
(with and without radiative cooling and star formation). Our 
results can be summarized as follows. 

• We show that gas distribution in simulated clusters has a 
rather spherical shape at large cluster-centric radii with 
the average axis ratios of ^ 0.8 at r > O.Srsoo. This 
implies that the standard assumption of spherical sym- 
metry in analyses of ICM using X-ray and Sunyaev- 
Zeldovich observations should be quite accurate. 

• We show that baryonic dissipation makes gas distribu- 
tion more spherical at 0.1 < r/r^Qo < 1, where the axis 
ratios are larger by ~ 0.2 on average in the cooling 
(CSF) runs compared to the non-radiative (NR) runs. 
At small radii r/rsoo < 0.1, the short-to-long axis ra- 
tios c/fl are lower in the CSF runs compared to the NR 
runs by as much as ^ -0.5, but the intermediate axis 
ratios b/a are reduced to a much smaller degree. The 
CSF gas distributions are oblate in their centers reflect- 
ing the presence of cool gas supported by rotation. 

• We present predictions for X-ray ellipticity profiles of 
intracluster gas based on mock Chandra X-ray maps 
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Figure 8. Ellipticity profiles averaged over mock X-ray maps of simulated clusters in the x-, y-, and z-projections (solid, dotted, and dashed lines, respectively). 
The blue lines represent the profiles for the CSF simulations, and the red lines for the NR runs. The left panel shows the profiles averaged over all clusters, 
the middle panel shows the profiles averaged over 21 projections, in which clusters appear relaxed, and the right panel shows ellipticity profiles averaged over 
projections in which clusters appear unrelaxed. The errorbars show the Ict errors of the mean values for the ellipticity profile for the x-projections. The errors are 
similar for the profiles in the y- and j-projections. 



of our simulated clusters. We show that the effect 
of cooling on the ellipticity of gas is consistent with 
the three-dimensional results. Specifically, the ellip- 
ticities at larger radii (0.1 < r/rsooe ^ 1) decrease for 
the CSF clusters compared to NR clusters reflecting 
the rounder gravitational potential of the CSF clusters. 
At r < O.lrsoo, ellipticities in the CSF clusters increase 
with decreasing radius. 

• NR clusters exhibit the opposite trends. NR clusters are 
more triaxial than the CSF clusters at r > O.lrsoo, but 
become considerably rounder at r < O.lrsoo- The latter 
trend is due to random gas motions, which are present 
at all radii but become considerable compared to the 
thermal energy of the gas (as reflected by increasing 
o-gas/v'circ ratio). 

• Our results indicate that observed ellipticity profiles of 
X-ray clusters can be used to constrain both the amount 
of cooling in the last several billion years of cluster evo- 
lution and presence of significant random gas motions. 

There are several issues to bear in mind when interpreting 
our results. First is the issue of overcooling in galaxy for- 
mation simulations. Our CSF simulations suffer from over- 
cooling in the cluster cores (r < O.lrjoo), such that the effect 
of halo contraction in response to the formation of the cen- 
tral galaxy is likely overestimated. The dissipational effect on 
gas shapes that we present can therefore be considered as an 
upper limit. At the same time, some amount of cooling must 
have occurred because galaxies do exist in clusters and the ob- 
served ICM gas fractions are signifi cantly below the expected 
values (e.g., Kravtsov et al.|[2009| . We therefore expect the 
ellipticity of the intracluster gas m real clusters to be between 
our results for CSF and NR runs. 

Although our cluster sample is not drawn to sample the 
mass function of clusters, the mass-dependence of elliptici- 
ties is very weak (Figure[4]i, at least in the range of masses we 
probe (3.5 x lO'^r'M© < Mgoo < 9 x lO'^/i-'M©) and our 
results should therefore not be biased significantly. 

A more significant concern is the possible bias due to the 
fact that our cluster sample was simulated assuming ug = 0.9, 



while the most recent estimates indicate erg = 0.80 ± 0.02 
(Vikhl inin et al. 1 2009; Jarosik et al. 2011). The effect of a 
lower cT g is that ha los would form later, leadin g to higher ellip- 
ticities ( [Allgood et al.^200 6). However, Maccio et al.|p008| l 
show that effect of changing o-g from 0.9 to 0.8 changes aver- 
age c/fl ratios of halos by only w 0.03, considerably smaller 
than the differences discussed in this paper. Moreover, both 
the CSF and NR clusters would be affected by the difference 
in cosmology and it is plausible that they would be affected to 
a similar degree. 

Another potential concern is that AGN heating is not in- 
cluded in our simulations though evidence of such heating 
is abundant i n real clusters, especially in r elaxed cool-core 
clusters (e.g., [McNamara & Nulsen 2007\ . Jets and bub- 
bles inflated by the AGNs can potentially change the shape 
of the gas significantly in the cluster core. We expect that the 
AGN would also change the gas shape in the cluster outskirts 
by suppressing the overall amount of cooling throughout the 
cluster formation. 

Despite these caveats, a few general implications can be 
drawn from our results. First, as dissipation makes the gas 
more spherical, systematics of observable quantities inte- 
grated along the line of sight such as Fsz, or any deprojected 
quantity that relies on the assumption of spherical symmetry, 
should be reduced compared to conclusions one could draw 
from dissipationless simulations. Another implication is that 
the shape of gas can be different from the shape of the po- 
tential in the cluster core r < O.lrsoo, and in the cluster out- 
skirts r > rsoo- The difference in shape between gas and po- 
tential can be attributed t o dev iation from hydrostatic equilib- 
rium as shown in Section|3.3| In the intermediate radial range 



0.1 < r/rsoQ < 1 , however, we have shown that the gas shape 
generally coincides with that of the potential, and therefore 
the shape of the potential can be inferred from the shape of 
gas. If there is an independent way of determining the shape 
of gravitational potential, e.g. by gravitational lensing, one 
may be able to constrain the amount of gas motions by com- 
paring the shape of the gravitational potential and the shape 
of gas. 

However, even without independent information about the 
potential, our results indicate that elUpticities derived from X- 
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ray images alone can constrain the amount of cooling and the 
presence of random gas motions. This is because these effects 
result in a rapid change of ellipticity with decreasing radius at 
r < O.lrjoo but of an opposite sign. We will present compar- 
isons of gas ellipticities in simulations and observations from 
Chandra and ROSAT in a separate companion paper (E. T. 
Lau et al. 2011, in preparation). 
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